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Abstract. Recent observations have shown that (some) metal line ab- 
sorption systems in QSO spectra arise from gas in galactic halos with 
radii much larger than the optical radii of galaxies. We show that the 
observed galaxy-absorber connections are natural results of galaxy for- 
mation in hierarchical cosmogonies, where galaxies form from gas cooling 
and condensation in dark matter halos. 



1. Introduction 

There is now much evidence that metal line absorption systems in QSO spectra 
arise from gas in galactic halos. Recent imaging and spectroscopic observations 
can identify directly, at low redshift {z <^ 1), the galaxies associated with the 
Mgll absorption systems (e.g. Bergeron, Cristiani, & Shaver 1992; Bechtold & 
Ellingson 1992; Steidel 1995; see Steidel 1998 for a recent review). The properties 
of such galaxy/ absorber systems can be briefly summarized as follows: (i) Most 
Mgll systems with rest-frame equivalent widths W > 0.3A are associated with 
galaxies; (ii) Most galaxies with (X-band) luminsity Lk ;> 0.05LJ^ have gaseous 

halos to produce Mgll systems with W > 0.3A; (iii) The typical absorption 
radius of a galaxy with L K ^ 0.05L* K is R a ~ 35 /i _1 Kpc(L^/L^) - 15 , much 
larger than the optical radius; (iv) Galaxies with Lk <J 0.05L* K have much 
smaller absorption cross sections than the brighter ones; (v) For bright galaxies, 
the absorbing clouds has roughly spherical distribution, with a covering factor 
^ 1 within R a , as indicated by the fact that few interlopers are observed; (vi) 
The absorbing galaxies show wide range of colors, from that of late type spirals 
to that of early type ellipticals, indicating that bright galaxies of all types possess 
similar gaseous halos; (vii) The observed absorption line systems show velocity 
structures, with typical velocity spread of ~ 50-200 km s -1 (e.g. Petitjean & 
Bergeron 1990; Churchill 1998), consistent with the absorption clouds moving 
in the potential wells of galaxies; (viii) The absorbing gas is enriched, as heavy 
elements are observed. 

Despite of the many observational facts about the relations between galax- 
ies and the metal line absorption systems, relatively few theoretical work has 
been done. In particular, the importance of such observational results to our 
understanding of galaxy formation has not been fully explored. 

In hierarchical clustering models of galaxy formation, it is assumed that 
galaxies form at the centers of dark matter halos as gas collapses, is shocked, 
and cools. Although our knowledge of how galaxies form is still incomplete at 



present time, what is certain is that at some time in the past gas had to dissipate 
and move from the outer parts to the inner parts of the dark matter halos around 
galaxies, and in this process it must have produced absorption line systems in 
any quasar line-of-sight intercepting galactic halos. It is clear that in order to 
understand galaxy formation we must first understand the physical state of the 
gas out of which galaxies form, and that absorption line systems are an excellent 
observational probe to any gas dissipating in halos. Thus, the observational 
evidence that the high column density absorption systems at z <J 1 indeed arise 
in gaseous halos around galaxies should be a landmark for our understanding 
of galaxy formation. Based on these considerations, we have proposed a simple 
two-phase model for the gaseous structure of galactic halos (see Mo, 1994; Mo 
& Miralda-Escude 1996). The model is based on current scenario of galaxy 
formation in hierarchical clustering models, and so allows us to understand the 
absorption line systems in the general framework of structure formation in the 
universe. It was found that the observed galaxy-absorber connections are natural 
results of galaxy formation in hierarchical cosmogonies. 

2. The Structure of Gaseous Galactic Halos 

During the process of galaxy formation, intergalactic gas collapse and move 
inward through the extended dark matter halos that are observed around present 
galaxies. A halo of hot gas at the virial temperature will form as the kinetic 
energy of the infalling material is thermalized in shocks. In low mass halos, the 
cooling time of this hot halo gas is short compared to the dynamical time if it 
contains all the accreted baryons; therefore, the gas must cool until the density 
of hot gas decreases to a level such that the cooling time is similar to the age 
of the system. Due to the increased cooling rate as temperatures drop, a two- 
phase medium will naturally form in these conditions of rapid cooling. Clouds 
photoionized by a radiation background will be maintained at an equilibrium 
temperature of about 10 4 Kelvin, in pressure equilibrium with the hotter halo 
gas. These clouds could form from inhomogeneities in the halo gas, or from 
the ram-pressure stripped interstellar medium of satellite galaxies. Due to their 
higher density, clouds will fall to the center of the halo where they may form 
globular clusters or settle into gaseous disks. In such a scenario, modeling the 
structure of gaseous galactic halos involves both the formation of dark matter 
halos and the state of gas in dark matter halos. The following is a sketch of the 
main ingredients involved in such a modeling. 

2.1. Dark halos 

The properties of dark matter halos in hierarchical clustering cosmogonies [such 
as cold dark matter (CDM) models] are relatively well understood from analytic 
models and N-body simulations. Dark matter halos can be described roughly as 
singular isothermal spheres, with density profiles given by 
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where V c is the circular velocity, M is the mass, and r v is the virial radius, of 
the halo. From spherical collapse model, the virial radius can be defined as 
r v = V c /WH(z), where H(z) is the Hubble constant at redshift z. 

For a given cosmogonic model, one can also estimate the mass function of 
dark matter halos, n(M; z)dM , which gives the comoving number density of 
dark matter halos, with mass in the range M — > M + dM, at a given redshift z. 
This function can be estimated either from N-body simulations, or from analytic 
models such as the Press-Schechter formalism. Thus, the number of halos which 
may host absorption line systems can be obtained. 

As demonstration, results in the following are shown only for the standard 
CDM cosmogony with Qq = 1, h = 0.5, and a% = 0.67. 

2.2. Gas in Dark Matter Halos 

Before the collapse of a dark matter halo, gas on a mass shell moves together with 
dark matter particles until the velocity of the mass shell reaches the sound speed 
of the gas interior to the mass shell, where the gas is shocked while dark matter 
particles continue to collapse and virialize. Ananlytic models (e.g. Bertchinger 
1989) and numerical simulations (e.g. Evrard 1990) show that, if cooling is not 
efficient then the gas in a dark halo will be shock-heated to the virial temperature 

T„ = M^, (2) 

with a profile similar to that of the dark halo. Because the gas density is higher 
near the center of a halo, radiative cooling of the gas is more efficient there. We 
define a cooling radius r c by 

fp ™ = 21^' (3) 

where A(T) is the cooling rate, fx is the average mass per particle, / is the 
gas mass fraction and tu ~ t/2 is the time between major mergers. Thus, gas 
located at r < r c can cool effectively before the halo merges into a larger system, 
while that at r > r c cannot cool and should retain its original state. Based on 
this consideration, we model the density profile of the hot component as 

PhotM = - J V ( C r, (4) 

47rGr(r + r c ) 

so that Phot( r ) f p{ r ) for r»r c and Phot( r ) oc 1/r for r r c . The tempera- 
ture of this phase is assumed to be at Thot ~ T v . These behaviours of />hot( r ) are 
similar to those obtained from the self-similarity solution given by Bertchinger 
(1989). When r c > r v , the cooling time at the virial radius is already short 
compared to tM- In this case we assume r c = r v . 

The accreted gas that does not stay in the hot phase and cools must fall 
through the hot halo in the form of photoionized clouds. We refer to the gas in 
these clouds as the "cold phase". The rate at which cold gas accumulates in a 
halo is determined by both gas infall and gas cooling, and is roughly 

fV c 2 

Mcoi d ~ — —r min ; r min = mm[r v ,r c ]. (5) 



We assume the mass flow rate to be M = M/tM, an d that the clouds move to 
the halo center with a constant velocity v. Assuming also spherical symmetry 
for the gas distribution, we can write the density of the cold g function 
of the distance r to the halo center as 

M . . 

Pcold(r) = (6) 

In the absence of a hot phase, the cold gas is in free fall and v must be of the 
order of the virial velocity V c . The friction of hot gas may cause cold clouds 
to move at a terminal velocity which is smaller than V c . However, since halos 
on galactic scales have cooling times that are similar to the dynamical times, 
the inflow velocity should not be much smaller than the virial velocity This is 
also consistent with the observed velocity structure in Mgll systems. The cold 
clouds are assumed to be at a temperature T c i ou d ~ 10 4 K, they are also assumed 
to be in pressure equilibrium with the hot gas, so that the density of a cloud 
(assumed to be uniform) at a radius r is 

Pdoud ~ Phot(r)T hot /T doud . (7) 

2.3. Cloud properties 

In the model outlined above, cold clouds are moving in diffuse substrates of hot 
gas. Various physical processes can act to affect the formation and destruction 
of the clouds. These processes constrain the masses and densities of the clouds. 
When the mass of a cloud is too large and its density too high, the cloud will 
become gravitationally unstable as it sinks towards the halo center. This hap- 
pens when the cloud mass exceeds the Jeans mass. Once a cloud is accelerated 
to a terminal velocity at which the drag force of the hot medium is comparable 
to the self-gravitating force on its surface, the cloud may become hydrodynam- 
ically unstable. Small clouds are also unstable against heat conduction. Based 
on these considerations, we found that the masses of cold clouds are confined to 
a narrow range, M c \ oud ~ 1O 6 M0. The typical size of such a cloud is ~ 1 Kpc. 
The covering factor of such clouds is ~ 1 within ~ 30 /i -1 Kpc for halos with 
V c ~ 200 km s _1 . 



3. Implication for QSO Absorption Line Systems 

To make connections between gaseous halos and QSO absorption line systems, 
we need to calculate the column densities of ions (of some species) at different 
impact parameters from a halo. We assume that the gas is ionized by the UV 
background radiation with a flux J(v) which we parameterize as 

J(v) = J-2i{z) x 10~ 21 (—) ©(^ergcm^sr^hz^sec- 1 , (8) 

where z^m is the hydrogen Lyman limit frequency, J-2i(z) gives the redshift 
dependence of J{v), and describes departures of the spectrum from a 

power-law. A typical spectrum can be represented approximately by taking 
a = 0.5, J-2i(z) = 0.5 for z > 2 and J_ 2 i(z) = 0.5 x [(1 + z)/3} 2 for z < 2, and 
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including a break in the spectrum at v± = 4 Ryd (due to continuum absorption 
by Hell), with @(v < 1/4) = 1 and 9 4 = Q(u > u A ) = 0.1. 

To consider metal line systems, we also need to assume a metallicity, Z, 
of the gas. The exact level of chemical enrichment in galactic halos is not 
clear. In the disk-spheroid model of our Galaxy, the metal output from stars in 
the spheroid component can enrich the disk gas to a level of about O.2Z (see 
Binney & Tremaine 1987, §9.2). In this model, the gas in the halo was initially 
enriched to such a level before it settles into the disk. In the galaxy formation 
model of White & Frenk (1991), gas in halos with V c ^> 100 km s^ 1 could also 
be enriched to a level of 0.1-0.3Z Q (see their Fig. 4), because of star formation 
in the progenitors and that in halos themselves. Unfortunately, the detail result 
depends crucially on how well the output metal is mixed with halo gas. In our 
model we treat Z as a free parameter. 

With this assumptions, the number density of any species of ions can be 
obtained as a function of the radius from the center of a halo with circular 
velocity V c at redshift z. As a result, for single halos we can obtain Nx(D), the 
column density of species X expected at an impact parameter D, and D(Nx), 
the impact parameter at which the column density is Nx- Given the number 
density of dark halos as a function of V c and z (e.g. from the Press-Schechter 
formalism), we can also obtain the statistical properties of absorption systems. 

The left panel of Figure 1 shows the HI column density Ahi versus impact 
parameters for halos at z = 0.5 with various V c . The value of z chosen here is 
to match the mean redshift of the Mgll systems for which associated galaxies 
are identified (e.g. Stcidel 1995). As one can see, Ahi decrease with D. This 
is because Ahi is proportional to the total column density of cold gas, which 
is higher for a smaller impact parameter, and to the pressure (since neutral 
fraction increases with pressure), which is higher near the center. The figure 
shows that the impact parameter for producing a Lyman limit system (LLS) 
with A H i £ 10 17 cuT 2 is about 40/i _1 Kpc at V c ~ 200 kms" 1 . In halos with 
V c ^ 250 km s -1 , this impact parameter is not much larger because the cooling 
radius is reached, and we assume that no clouds are formed outside r c . 

The shape of the Mgll column density profile (shown in the right panel of 
Fig.l) is similar to that of HI, because A^gii/Ani is roughly a constant for the 
photoionization model considered here. The observed samples of Mgll absorp- 
tion lines select systems where the Mgll line is above a threshold of equivalent 
width. We present results for the number of absorbers and the impact parameter 
distribution assuming a fixed threshold AM g n = 2 x 10 13 cm~ 2 (reasonable for 
W = 0.3A), although we should bear in mind that the average cloud covering 
factor can also affect the observed equivalent widths. The left panel of Figure 2 
shows the results for the impact parameter in halos of different V c , at redshift 
z = 0.5. The UV flux spectrum used here has parameters J_2i = 0.1, a = 0.5 
and Q(v > 1/4) = 0.1. At V c ^ 150 kms -1 , D decreases rapidly with decreasing 
V c , due to the faster cooling rate and the lower pressures that are implied. The 
typical impact parameter is about 30/i -1 Kpc at V c ~ 200 km s -1 . In halos with 
V c ;> 250 km s -1 , this impact parameter is not much larger again because the 
cooling radius is reached and no clouds are assumed to form outside the cooling 
radius. This result of impact parameter as a function of V c is similar to the 
observational result of Steidel et al. (see Steidel 1995). The observation shows 
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Figure 1. Left panel: HI column density as a function of impact 
parameter D, for halos at z = 0.5 and with various V c . The model 
parameters are / = 0.05, Z = 0.3Z© and v = V c . The horizontal line 
indicates Nm = 10 17 cm~ 2 , above which an LLS is produced. Right 
panel: Mgll column density as a function D. Results are shown for 
two models of UV flux, one is a power law with power index a = 0.5, 
the other has the same power index but with a break at 4 Ryd so that 
64 = 9(f > v*) = 0.02. In both cases J_ 2 i = 0.1. 




Figure 2. Left panel: The impact parameter D, at which Nyiga = 
2 x 10 13 cm~ 2 , as a function of halo circular velocity. Results are show 
for halos at z = 0.5. The dotted line shows the relation D oc V c 6 . 
Right panel: Cumulative cross section of Mgll systems at z = 0.5, 
with iV~Mgii > 2 x 10 13 cm 2 , given by halos with circular velocity V c > 
Vq, as a function of Vq. 
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Figure 3. Left panel: Number of absorption systems per unit 
redshift, for HI systems with JVhi > 10 17 cm~ 2 , Mgll systems with 
A^Mgii > 2 x 10 13 cm~ 2 , and CIV systems with iV C rv > lutein" 2 . Solid 
curves assume / = 0.05 and v = V C with Z = 0.3Z©; dashed curves as- 
sume the same / and v but with Z = 0.1Z & . Right panel: The impact 
parameters at which AT MgII = 2 x 10 13 cm -2 (solid), iVciv = 10 14 cm -2 
(short dashed), and Afovi = 10 14 cm~ 2 (long dashed), as a function of 
V c for halos at z = 0.5. 

that the maximum impact parameter D of their Mgll systems to the absorbing 
galaxies (having typically Lk ^ 0.05-L^) increases slowly with the K-band lu- 
minosity Lk as D oc L° K 15 . Assuming the Tully-Fisher relation Lk oc V c a , we 
have D cx V c °' & , which is shown in the left panel of Fig. 2 by the dotted line. We 
see that such a weak increase of D with V c can be accommodated in our model, 
arising from the value of the cooling radius and the derived gas pressures in dif- 
ferent halos. The very small impact parameters in low V c halos agrees with the 
fact that Mgll absorption systems are not commonly found in galaxies fainter 
than about 0.1L* (Steidel 1995). 

To see how different halos contribute to the Mgll absorption, we plot in the 
right panel of Fig.2 the cumulative number of Mgll systems per unit redshift 
at z = 0.5, produced in halos with V c < Vq. The figure shows that most Mgll 
systems are produced in halos with V c = 150-300 km s _1 . The median V c is 
about 200kms~ 1 . The contribution to the total cross section made by halos 
with V c <; 150 km s" 1 is small, which again suggests that Mgll systems should 
not be commonly found in the halos of dwarf galaxies. 

The left panel in Fig.3 shows the total number of Mgll systems and the 
evolution with redshift, compared to the number of LLSs, for two values of the 
metallicity. The relative number of Mgll systems to LLSs is correctly predicted 
when the metallicity is between Z = 0.1-0.3, (compared with the observational 
results in Steidel & Sargent 1992), and the evolution is also as observed if the 
metallicity in these halo clouds does not decrease very fast with redshift. For 
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comparison, the total number of CIV systems in our model is also shown in the 
panel; the ratio to the number of Mgll systems, and the evolution with redshift, 
is similar to what is observed. Fig. 3 (right panel) also compares the impact 
parameters of CIV systems with those of Mgll systems in different halos. In 
small halos, D is much larger for CIV systems than for Mgll systems. In fact, 
the cross-section weighted average of V c is only about 100 km s -1 for the CIV 
systems, but as large as about 200kms _1 for Mgll systems, in the standard 
CDM model considered here. Our model thus predicts that, while Mgll systems 
are mostly associated with bright galaxies, many CIV systems should be found 
to be associated with small galaxies where A^Mgii is low. For the same reason, 
OVI systems should neither be commonly found as photoionized clouds in halos 
with V c ^ 150 km s -1 . Some of such highly-ionized systems may arise from colli- 
sionally ionized gas (e.g. in the hot phase) or from an intermediate temperature 
phase between the hot medium and the cold clouds. 

4. Discussion 

We have demonstrated that the simple halo model discussed in §2 can indeed re- 
produce the main properties of the observed galaxy-absorber connections. This 
is important, because it means that the observed absorption systems can be un- 
derstood in the current framework of galaxy formation. Much theoretical work 
remains to be done to see the plausibility of our model and to make it more 
predictive, mainly to understand the physical processes determining the rate of 
formation and destruction of gas clouds moving through the halo, the distri- 
bution of cloud masses and velocities, etc. Further observations of absorption 
line systems and their relation to galaxies should open a new era in our under- 
standing of how galaxies form, as the physical conditions of the gas in the halos 
around galaxies, and the dependence of these conditions on galaxy luminosity 
and morphology are unravelled. 
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